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Abstract 

Although the quantum bound-state energies may be generated by the so called 
PT— symmetric Hamiltonians H = VH'^V 7^ if^ where V is, typically, parity, the 
spectrum only remains real and observable (i.e., in the language of physics, the 
"PT— symmetry remains unbroken) inside a domain V of couplings. We show that 
the boundary &D (i.e., certain stability and observability horizon formed by the 
Kato's exceptional points) remains algebraic (i.e., we determine it by closed formu- 
lae) for a certain toy-model family of A^— dimensional anharmonic-oscillator-related 
matrix Hamiltonians H^^^ with A^ = 2, 3, . . . , 11. 



1 Introduction 



According to the abstract principles of Quantum Mechanics, the observable quan- 
tities (say, the spectra of energies Eq < Ei < . . . of bound states) should be con- 
structed as eigenvalues of a certain self-adjoint operator H = H'^ acting in some 
physical Hilbert space of states Ti. Fortunately, the full and impressive generality 
of this formulation of the theory is rarely needed in its concrete applications. Most 
often, the Hilbert space is being chosen in its most common representation Ej2{JR) 
with elements \ip) representing the square- integrable complex functions of a single 
variable x interpreted as a coordinate of a (quasi)particle. 

The most common version of the Hamiltonian H composed of its kinetic and 
potential-energy parts leads to the constructions of the energies via a suitable phe- 
nomenological potential V{x) entering the ordinary differential Schrodinger equation 

-^^i^{x) + V{x)ij{x) = Ei;{x). (1) 

In 1998, Bender and Boettcher [Ij demonstrated that the spectra {En} can remain 
real, amazingly, even if the potentials themselves become complex. This attracted the 
attention of physicists of different professional orientations ranging from supersym- 
metry |2j and field theory [3j to cosmology [4J and even magnetohydrodynamics [3]. 

Although the Bender's and Boettcher's observation has been supported by many 
subsequent studies and concrete examples [6l [TJ [8] it still looked like a paradox as 
it immediately implied that H ^ W in Ej2{M). Fortunately, the resolution of the 
paradox proved rather easy [9j. It was sufficient to imagine that the Hamiltonian H 
can remain tractable as self-adjoint in another Hilbert space of states 7i ^ Ej2{M). 
A closely related historical paradox is that the usefulness of transfer of H ^ W to 
another Hilbert space Ti (where it becomes self-adjoint) has already been known and 
even applied successfully in nuclear physics in 1992 [TU] . Still, it took time to clarify 
this parallelism [11] (cf. also the comprehensive review paper [12] for further details). 

On the level of the elementary ordinary differential example ([1]), one of the most 
important formal observations has been made in ref. [B]- It's authors Dorey, Dunning 
and Tateo noticed and described a deep nontriviality of the problem of the reality 
of the spectrum {En}- After they added some other free parameters (the physical 
meaning of which is not too relevant for our forthcoming argumentation) and after 
they replaced the Bender's and Boettcher's if = _|_ ^2 (^^^y with 5 > by a 
two-parametric family of their generalized Hamiltonians H^^^^\a,£), they revealed 
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that the reahty of the spectrum {En{a,i)} only takes place inside a certain domain 
of parameters [let's denote it as T^la, i)] with a highly nontrivial, spiked shape of its 
boundary dV{a,i). 

We found the latter observation extremely challenging, important, interesting 
and inspiring. The very existence of a finite boundary of the domain V{a,i) ^ 
represents one of the main differences of the model from the differential Schrodinger 
eqs. ([1]) which are self-adjoint in the standard Hilbert space Ki2{JR) and which require 
that the potential V{x) remains real. In the new framework, the change of space 
Ej2{M) — * H implies that only the new potentials V{x) with parameters inside V^a, i) 
may consistently be chosen as complex |12j . 

The main motivation of the detailed study of the domains exemplified hj V{a, £) 
lies in their obvious practical relevance. Virtually no problems emerge in the current 
textbook scenario where the explicit specification of the physical domain of parame- 
ters V is usually trivial. The standard choice of the space Ej2{M) and of a manifestly 
self-adjoint differential Hamiltonian H = (depending on some J real couplings or 
other free parameters) usually enables us to work with the elementary, unrestricted 
domain T) = of these parameters. In contrast, the necessary determination of 
dV (which plays the role of a certain horizon of the observability and of the stability 
of the quantum system) proved fairly difficult even in the comparatively elementary 
example H^^^^){a,£) ^ [H^^^^\a,£)f as studied in ref. [6]. 

In the widespread terminology coined in refs. [H [9] and employed also in the 
review [12], the set V{a,i) of parameters could be called the domain of the so called 
unbroken PT— symmetry of the system in question The scope of this specification is 
adapted to the models ([1]) - usually, one stays inside Ej2{M) and specifies V as the 
parity reversal and T as the time reversal. Then one requires the PT— symmetry 
of the Hamiltonians H which means that VTH = HVT PQ [T3] . The choice of the 
parameters inside V (i.e., the physical requirement of the measurability and reality 
of the spectrum) is then, finally, rephrased as the so called PT— symmetry of the 
wave functions (the review [12] can and should be consulted for all these details). 

The nontriviality and the non-smooth, spiked shape of the curve &D{a,i) as 
found in [S] opens the question of a generic characterization of the geometry of 
dV in the less specific setting. In what follows, we intend to review, extend and 
complete our results in this direction as published in the papers jH] - [21]. In these 
papers we tried to bridge the apparent gap between the methods aimed at differential 
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Schrodinger operators and at their alternative matrix representations with H ^ W . 
A supplementary though still sufficiently appealing physical background of our series 
of studies has been found in the anharmonic-oscillator problem with the specific and 
popular H = H^^^'~'^ ^ H'^ where V{x) = V^^^'~'\x) = ^\gx^ or where Vix) has 
been generalized to an arbitrary real polynomial in the purely imaginary variable \x. 

The results of refs. [H] - [2T] (with particular attention paid to the constructions 
of the domains "D) will be reviewed here in sections |2] and [31 Their extension and 
completion will be described in sections H] and [51 Section [6] is a summary. 

2 Matrix models with small dimensions 

In the review [lOj written in the context of the so called interacting boson models in 
nuclear physics, Scholtz et al emphasized that once a non-unitary Dyson's fermion- 
to-boson mapping is used, the calculation of the energies gets simplified while the 
price to be paid still proves reasonable for matrices. Indeed, in order to restore the 
Hermiticity of the Hamiltonians, just the finite- dimensional physical Hilbert space 
7i had to be reconstructed ad hoc. This observation of merits of working with 
finite-dimensional matrices also formed a key encouragement of our forthcoming 
considerations. 

2.1 Variational inspiration: truncated Hamiltonians 

By assumption, our starting point given by Schrodinger eq. ([T]) with H = if (^^'^) 
remains simple and tractable in the usual Hilbert space Ej2{M). At the same time, the 
necessary transition to the physical Hilbert space Ti (where our complex anharmonic- 
oscillator becomes self-adjoint) may remain complicated (interested readers 
can find all details in the literature cited in [I2|)- In such a situation, the use of the 
basis {plus its subsequent variational truncation) has been recommended in all the 
series of our studies Jl] - |2l]. In this way we were able to reduce the differential- 
operator Hamiltonians H^"^^^^ in eq. ([1]) to the sequence of their partitioned matrix 
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approximate forms 



A 


B 


C 


D 











Hk,i 


Hk,k 


Hk,k+i 


Hk,k+k' 


Hk+1,1 


Hk+i,k 


Hk+i,k+i 


Hk+i,k+k' 


Hk+k',1 ■ 


■ Hk+k',k 


Hk+k',k+i ■ 


■ Hk+k',k+k' , 



The truncated basis {\n, ±)} has been numbered by the excitation quantum numbers 
n = 0, 1, . . . and it has been partitioned with respect to the parity ±. For our special 
anharmonic cubic H^^^'~'^ this imphed that A and D were real and diagonal, with 
^k,k = 4:k + l and D^ i^ = 4/c + 3 in suitable units. In parallel, all the necessary matrix 
elements of the purely imaginary matrices B = iB' = C* = i(C")* (where * means 
transposition) have to be evaluated numerically. This is partially simplified by the 
observation that we may restrict our attention to the purely real matrices since 
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C 


D 




A 


B' 


-a 


D 




Unfortunately, our attempts to construct the domains V in closed form were only 
successful at N = K + K' = 2 [14j, atN = K + K' = 3^^ and, with certain 
surmountable difficulties, at N = K + K' = 4 [17] . 



2.2 Tridiagonal chain models in the strong-coupUng regime 

The latter results helped us to imagine that the variations of the vast majority of the 
(real) matrix elements of B' = {C'Y did not lead to any really significant changes in 
the overall structure of the respective boundaries d'D^^\ For this reason we further 
re-arranged the basis and reduced the variability of the submatrices B' and C by 
setting many of their matrix elements equal to zero. As a net result of all these more 
or less natural simplifications we arrived, in ref. [T5j, at another family of the matrix 
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toy models exhibiting the general tridiagonal structure 



-(iV-1) 

-gi -{N - 3) 

-92 



92 



N -5 
-92 



92 

N -3 
-9i 



91 



(2) 



At the dimensions = 2 J or = 2 J + 1 these models depend just on a J— plet 
of real couplings gi,g2, ■ ■ ■ ,gj- One of the most important formal merits of these 
models is that at any dimension A^, all the domain X)(^) lies inside a finite hypercube 
[ig. The boundary dV^^^ itself can be characterized by its strong-coupling maxima 
which were obtained in the following closed form, 



n 



1,2, 



J. 



(3) 



Although the strong-coupling result ^ looks easy, its derivation required extensive 
computer-assisted symbolic manipulations. Via a nontrivial extrapolation guess- 
work we revealed that geometrically, the horizons dV^'^^ are (hyper) surfaces with 
protruded spikes called extreme exceptional points, EEPs. This intuitive picture has 
been complemented by the more quantitative descriptions of dV^^^ in [20] . It 
was based on the strong-coupling perturbation ansatz using an auxiliary, formally 
redundant small parameter t, 



9n = 9 



(4) 



This ansatz extrapolates the rigorous J < 2 fine-tuning rules as derived in refs. [14] 
m to all J. 



2.3 Secular equations 

Once we choose N = 2J or N 



2J + 1, abbreviate E'^ 



s and, at all the odd 
0, we find that 

all the secular equations det (^H'^^^ — = have the same polynomial form. 



dimensions A^ = 2 J + 1, ignore the persistent energy level 



J 




,J-2 



Q 




,J-3 



R + 



0, 



(5) 
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At all J and A^, the coefficients P,Q,R,... should be understood as real polynomial 
functions of the J— plets of squares gl of our real matrix elements. Once all the ener- 
gies are assumed real (i.e., equivalently, once all the roots Sk of eq. ([5]) happen to be 
non-negative), we immediately deduce the following relations tractable as necessary 
conditions imposed upon our coefficients in ([5]), 

J 



^ . Si + S2 + . . . + Sj > , 

Q = S1S2 + S1S3 + . . . + SiSj + S2S3 + S2S4 + . . . + > . 



(6) 




In the opposite direction, the set of the necessary inequalities P > 0, Q > 0, . . . is 
incomplete as it does not provide the desirable sufficient condition of observability. It 
admits complex roots s in general (take a sample secular polynomial (s^ + l)(s— 2) for 
illustration). This shows that our problem of the determination of the physical do- 
mains T'(^) of couplings is mathematically nontrivial even at the smallest dimensions 
N and J. 

For a given prototype Hamiltonian H^^^ and under the constraints the de- 
termination of the domain T'(^) = V (^H^^^^ is equivalent to the guarantee of the 
non-negativity of all the J roots Sk of eq. (Q. Keeping this idea in mind, the explicit 
forms of the corresponding sufficient conditions are to be given here for the ffist ten 
smallest matrix dimensions iV = 2,3,...,ll. 

3 The domains and p(2J+i). ^ brief review of 

the known results 

3.1 Methodical inspiration: the non-negativity of the root 
of eq. dSD at J = 1 

The ffist nontrivial illustration of the current Hermitian Schrodinger's bound-state 
problem is provided by the two-by-two real-matrix model 

H\il)) = E\il)) , H = H{a,b,d) = \ ^ ^ \ = W{a,b,d) 

b d 
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Its three-parametric spectrum is always real and, therefore, observable. 



E = E±{a,b,d) = - a + d ± ^ {a - dy + Ab^ 



For H = H{a, b, d) = W the three-dimensional physical domain T>{a, b, d) of param- 
eters giving real spectra coincides with all IB?. 

The parallel PT— symmetric two-by-two example is very similar. 



For each individual choice of the parameters a, b and d, the reality of the spectrum 
E'^{a., b, d) of the primed Hamiltonian H'{a, b, d) is fragile and it must be guaranteed 
and proved at a given triplet of parameters. The reality and stability of the primed 
system can only be achieved inside a perceivably smaller domain T>'{a, b, d) with the 
easily specified EP horizon. 



The interior of the non-compact manifold I'' (a, b, d) is specified by the single ele- 
mentary constraint 6g (— |a — (i|,|a — This may be interpreted as a fact that 
the variability of the parameters a and d is entirely redundant for qualitative con- 
siderations. It makes sense to get rid of them by the multiplicative re-scaling of all 
the parameters and by the subsequent shift of the energy scale leading to the generic 
choice of a = — 1 and d = 1. 

In the context of section 12. 3^ the latter reduction leads us to the linear version 
s — P = of secular eq. ([5]) at J = 1 which has the single root sq = P- The non- 
negativity of this root is equivalent to the non-negativity of the coefficient P. In 
terms of the single coupling gi = a available at J = 1, the necessary and sufficient 
criteria of the observability of H'^'^^ = H^'^\a) or if^^^ = H^^\a) read P^'^\a) = 
1 — > and P^^\a) = 4 — 2a^ > 0, respectively. In a way transferable to any 
dimension, the explicit definitions V'^'^^a) = (—1,1) and T)^^\a) = (— v^, v^) may 
be re-read as definitions of the corresponding EP horizons dT>^'^\a) = {—1, 1} and 



In the language of phenomenology, one notices an important complementarity 
between the parameter-dependence of the toy spectra E±{a,b,d) and E'^{a,b,d) at 
N = 2. In the former example, all of the energies E±{a, b, d) remain safely real. The 





dV^^\a) = {-V2, V2}. 
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second, primed model is less easy to deal with. There exists the whole set of the 
eligible two-by-two metric operators O = > which define the inner product in 
the corresponding two-dimensional toy Hilbert space Ti' (cf. [15j). Thus, the operator 
H'{a, b, d) represents an observable, in Ti' , in spite of its manifest non-Hermiticity 
in the auxiliary two-dimensional Hilbert space Ti (where the metric is the Dirac's 
simplest identity operator). With parameters inside V'{a,b,d) the primed model 
remains safely compatible with all the postulates of Quantum Mechanics, therefore. 

3.2 The non-negativity of all the roots of eq. ([5]) at J = 2 

At J = 2 the quadratic version s'^ — 2Ps + Q = 0of secular eq. ([5]) has two roots 
s-t = P ± — Q. These two roots remain real if and only if B = — Q > 0. 
In the subdomain of parameters where P > they remain both non-negative if and 
only if P > and Q > 0. We can summarize that the required sufficient criterion 
reads 

P > , P2 > Q > . (7) 

In an alternative approach, without an explicit reference to the available formula for 
s±, let us contemplate the parabolic curve y{s) = — 2 P s which remains safely 
positive, in the light of our assumption ([6]), at all the negative s < 0. This curve can 
only intersect the horizontal line z{s) = —Q at some non-negative points s > 0. 

The proof of non-negativity of all the roots of our secular equation degenerates 
to the proof that there exist two real points of intersection of the J = 2 parabola 
y{s) with the horizontal line z{s) (which lies below zero) at some s > 0. Towards 
this end we consider the minimum of the curve y{s) which lies at the point Sq such 
that ?/'(so) = 0, i.e., at Sq = P. This minimum must lie below (or, at worst, at) 
the horizontal line of z{s) = —Q < 0. But the minimum value of ?/(so) is known, 
y{P) = —P^. Thus, the condition of intersection y{so) < z{so) gives the formula 
P2 > Q. QED. 

It is amusing to notice that once eq. dHD holds, the inequality P^ — Q > is 
equivalent to the reality of the roots simply because P^ — Q = (si — S2)^/4. Even for 
some other two-parametric matrices, precisely this type of requirement is responsible 
for an important part of the EP boundary &D (cf. refs. [T6t [T7] for details). 
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4 The domains and I^^^J+i). Yiew results 

The determination of the physical horizons dV^^^ of our models H^^^ becomes a 
more or less purely numerical task at the very large matrix dimensions N [20] • In 
an opposite extreme, as we already noticed, the non-numerical exceptions have been 
found at = 2 [14J, at = 3 [16] and at the next two dimensions = 4 and 
N = b [T71 [TS] . Now we intend to complement these observations by showing that 
the closed-form constructions of the prototype horizons dV^^^ remain feasible up to 
the dimension as high as = 11. 



4.1 The non- negativity of all the roots of eq. ([5]) at J = 3 



Neither at A^ = 6 nor at A^ = 7 the sufficient condition of non-negativity of all the 
energy roots s is provided by the three necessary rules -P > 0, Q > and i? > of 
eq. Let us return, therefore, to the second method used in paragraph 13.21 and 
derive another inequality needed as a guarantee of the reality of the energies. In the 
first step one notices that all the three components of the polynomial 

y{s) = -SPs'^ + 3Qs = R, J = 3 

remain safely non-positive at s < 0. Whenever the roots are guaranteed real, their 
non-negativity s„ > with n = 1,2,3 is already a consequence of the three con- 
straints ([6]). A guarantee of their reality is less trivial but it still can be deduced 
from the shape of the function y{s) on the half-axis s > 0, i.e., from the existence 
and properties of a real maximum of y{s) (at s = s_) and of its subsequent minimum 
(at s = At both these points the derivative y'{s) = 3s^ — 6Ps + 3Q vanishes 
so that both the roots s± = P ± \/P'^ — Q of y'{s) must be real and non-negative. 
This condition is always satisfied for the real roots Sk of y{s) since 



54 L 



> 0. 



In the next step, the sufficient condition of the reality of the roots Sk will be un- 
derstood as equivalent to the doublet of the inequalities ?/(s_) > R and y{s^) < R. 
Here we may insert = 2Ps± — Q and get the two inequalities which are more 
explicit, 

2(p2_Q)s_ < PQ-R < 2{P'^-Q)s+. (8) 
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They restrict the range of a new symmetric function of the roots, 
PQ - R = ^ [siS2 (si + S2 - 2 S3) + S2S3 (s2 + S3 - 2 Si) + S3S1 (S3 + si - 2 S2)] . 
After another insertion of the known s± we arrive at a particularly compact formula 
2(p2 _ Q)3/2 > R-3PQ + 2P^ > -2 _ g)3/2 

or, equivalently, 

4 [p^-qY > [R-3PQ + 2P'^y . 

Due to the numerous cancellations the latter relation further degenerates to the most 
compact missing necessary condition 

3P^Q^ + 6RPQ>4Q^ + R'^ + 4RP^ . (9) 

Our task is completed. In combination with eqs. (jH]), equation plays the role of 
the guarantee of the reality of the energy spectrum. 

4.2 The non-negativity of all the roots of eq. ([5]) at J = 4 

In a search for the non-negative roots of the quartic secular equation 

det (i/^^'^) - E /) = X* -APx^ + QQx'^ -4Rx + S = y{x) + S = (10) 

we note that all the four A^— dependent coefficients P, Q, R and S evaluate as 
polynomials in the squares of the four coupling parameters Qk, k = 1,2,3,4. Once 
all these four expressions are kept non-negative, the curves y{x) and z{x) = —S do 
not intersect at x < 0. At x > they do intersect four times at x > (as required), 
provided only that the three extremes of y{x) can be found at the three non-negative 
real roots Xi^2,3 of the extremes-determining equation 

l/'(a;i,2,3) =4(x3 2,3-3Px?,2,3 + 3Qxi,2,3-^) =0. (11) 

In an ordering < Xi < X2 < X3 of these roots we arrive at the three sufficient 
conditions 

y{x,)<-S, y{x2)>-S, y{x3)<-S (12) 
guaranteeing that the parameters lie inside V^^^ or T)^^\ 
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All the three quantities Xk satisfy the cubic equation y'{x) = so that its pre- 
multiplication by x and enables us to single out the fourth powers of the roots, 

^1,2,3 = 3 Pxf 2,3 ~ 3 Q 3^1,2,3 + -^2:1,2,3 • 

In the other words, we can eliminate all the fourth powers of these roots from 7/(0:1,2,3) • 
This simplification reduces all the three items in eq. (|T2l) to the three polynomial 
inequalities of the third degree, 

-Pxfg + 3Qxl^ - 3i?si,3 + S <0, 

-P xl + 3Q xl - 3 Rx2 + S > 0. 

Repeating the same elimination of the maximal powers once more, we may insert 
^1,2,3 = 3 Pxf 2,3 — 3 Q Xi,2,3+-R and arrive at another equivalent triplet of inequalities 

-Bxj + 2B^/^Cxi<B^D, (13) 
-Bxl + 2B^/^Cx2> D, (14) 
-Bxl + 2B^^^Cx3< B^ D. (15) 

The old abbreviations S = _ g and 2 5^/2 (j ^pQ_ji pj^g ^ new one, 3 B'^ D = 
P R — S enable us to define ^1,2,3 := 2^1,2,3, /V^- This yields our final triplet of the 
very transparent quadratic-equation conditions 

- 2CY1 + D >0, (16) 

- 2CY2 + D <0, (17) 
Y^ -2CY3 + D>0. (18) 

The auxiliary roots Y± = C ± \JC'^ — D must be real and non-negative so that we 
must guarantee that D > and > D. The conclusion is that eqs. (fT6|l - (fT8|) 
degenerate to the four elementary requirements 

Yi<Y^<Y2<Y+<Y;. (19) 

Together with the inequalities P>0,(5>0and-l<C - + Q/B < 1 they 
form the final and complete algebraic definition of the domains V^^^ and V^'^^ . 

We can summarize that at J = 4, the feasibility of the non-numerical construction 
of the domains "D'-^-* and T'^^-' [determined by eqs. (jH]) and flTIJl) ] is based on the non- 
numerical solvability of the third-order polynomial equation (ITTl) . 
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4.3 The non-negativity of all the roots of eq. ([5]) at J = 5 

Let us finally proceed to H^^^ with N = 10 and/or = 11 which leads to the secular 
equations of the fifth degree, 



-5Px^ + 10Qx^ -lORx^ + 5Sx -T = y{x)-T = 0. 



(20) 



From our present point of view the problem of the construction of the respective 
horizons dV^^^ remains solvable exactly since the derivative y'{x) is still a polynomial 
of the mere fourth degree, 



The exact, real and non-negative values xi < X2 < x^ < x^ of the four roots of y'{x) 
(which determine the extremes of the function y{x)) may still be considered available 
in closed form. 

In a way which parallels our preceding considerations we may assume that the five 
A^— dependent non-negative coefficients P>0,Q>0,R>0,S>0 and T > obey 
also all the additional inequalities derived in the preceding sections. We may then 
treat our secular problem (!20!) as a search for the graphical intersections between the 
(nonnegative) constant curve z{x) = T and the graph of the polynomial y{x) of the 
fifth degree (which can only be nonnegative at x > 0). 

Inside the domain T)^'^\ the quintuplet of the (unknown but real and nonnegative) 
physical energy roots Xa, Xb, Xc, Xd and Xe has to obey the obvious intertwining rule 



The way towards the sufficient condition of the existence of the real energy spectrum 
remains the same as above, requiring 




(21) 



< Xa < Xi < Xb < X2 < Xc < Xs < Xd < X4 < Xe ■ 



y{xi) > T 



y{x2) < T 



y{x3)>T, y{xi)<T. 



(22) 



The lowering of the degree should reduce eq. fl22|) to the quadruplet 



w{Yi) < 



w{Y2) > 



w{Ys) < 



w{Yi) > 



(23) 



where the re-scaling a;i^2,3,4 = ^1,2,3,4 VB applies to the arguments of the brand new 
auxiliary polynomial function of the third degree in Y, 



w(Y) = Y^ -3CY^ + 3DY -G. 
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Besides the same abbreviations as above, we introduced here a new one, for PS—T = 
4 5^/^ G. The new and specific problem now arises in connection with the necessity 
of finding the three auxihary and, of course, real and non-negative roots of the 
cubic polynomial w{Y). Once we mark them, in the ascending order, by the Greek- 
alphabet subscripts, we should either postulate our (in principle, explicit) knowledge 
of their real and nonnegative values <Yp <Y^ or, in another perspective, we may 
recollect simply the above-derived conditions which restrict the range of the three 
coefficients C, D and G in the cubic polynomial w{Y). 

We may immediately conclude that the last feasible specification of the domains 
25(10) g^j^^ will be given by the following set of the inequalities, 

Yi<Y^<Y^<Yp<Y^<Y,<Y,. (24) 

This is the desired set of the missing algebraic formulae which complete the sufficient 
condition of the reality of the spectra. We can emphasize that at J = 5, the feasibility 
of the present non-numerical constructions of the most complicated five- dimensional 
though still algebraic domains V^^^ is related again to the most complicated though 
still non-numerical solvability of the extreme-determining fourth-order polynomial 
equation (12T|) . 



The series of the solvable models is, obviously, exhausted. Any attempted ex- 
tension of the recipe beyond = 11 would suffer from the necessity of using mere 
numerical auxiliary functions of couplings Qk representing the roots of the extreme- 
determining higher-order polynomials. 



5 The wedges of the hypersurfaces dV^^^ 

The existence of the algebraic formulae which determine all the boundaries dV^^^ 
up to = 11 opens a way towards a verification of the strong-coupling perturbation 
results of refs. fTEl and [Ldj. Hopefully, some other similar qualitative or geometric 
features of the observability horizons &D assigned to a given PT— symmetric Hamil- 
tonian H will be also detected or discovered via the reduction of the model to a series 
of its A^ by A^ approximations of the prototype form H^^\ 

In an alternative setting, we can fix the dimension A^ and search for the generic 
features represented by the given model H^^\ Moving beyond the perturbative 
framework, for example, our ansatz (jlj) could be then re-interpreted as a precise 
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change of the variables in the space of our free dynamical parameters. The redundant 
measure of smallness t may be now fixed arbitrarily. Instead of the couplings one 
can also decide to work with the free parameters 7^ or even with G^s. For illustration 
one may recollect the two-by-two Hamiltonian 



Z_ . ae(0,l) (25) 




with the two-point spectrum Ej^ = ±y/a . The variability of the new parameter 
coincides with T)^'^\a) = (0, 1) since there are no additional constraints. 

At the higher dimensions N, our horizons dV^'^^ may be shown to exhibit a 
generic hedge-hog-like shape as well as certain reflection symmetries. They allow us 
to restrict our attention to the subdomains of T)^^^ with the positive gkS, i.e., with 
the real quantities 'jj^^ = a, 7j^i = /3, . . . which should all remain non-negative and 
smaller than one, j^^^ G (0, 1). 

5.1 New forms of approximations 

Let us now return to the description of the structure of the boundaries of the domains 
p(^) at = 6 and = 7 by means of our key inequality ([8]). In this form of the 
rigorous guarantee of the reality of the energies at J = 3 we may set = P'^{B,Q) = 
B + Q, postulate B > 0, Q > and insert 



s± = ^B + Q±\^>0 

m eq. dHD. With an abbreviation C := (PQ - R)/{2B^/'^) > this gives the pair of 
inequalities S- < C \fB < s+ or, equivalently, 

TTq-1 < C < Jl + q + 1, g=^e(0,cx)). (26) 
* B 

Once we notice that PQ = Q y/B + Q = qB^I'^^X -|- q we may return from the 

auxiliary C to the original R and rewrite eq. fl2B]) as our final, perceivably simplified 

one-parametric constraint imposed upon the allowed range of the variability of the 

value of the polynomial i?, 

-(f-)v^.^.{L„_j;:'. (^^) 



In such a reparametrization of the physical domains D^^''^) we use S > and Q > 
and employ our final form of the two-sided inequality ( 1271) as a definition of the 
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allowed range of the remaining quantity R. The latter inequality if particularly 
strong at the smallest ratios q = Q/B. As long as 

V2 Jy 8^ 8^ 128^ 64^ 1024^ v"^ / ' 

we see that the smallness of q implies the second- order smallness of R. This means 
that in the regime where Q <^ B, the three-dimensional physical domains are 
very narrow in their third dimension represented by R. 

5.2 The pairwise confluences of the levels at J = 3 

For the sake of definiteness let us choose just = 6 and abbreviate 

gi = c < -\/5 , 92 = b < 2 V2 , 93 = a < 3 

in the six-by-six version of matrix ([2]). It is easy to deduce that the domain T)^^^ is 
circumscribed by the ellipsoidal surface given by the equation 

P = - (a^ + 2 6^ + 2 - 35) /3 = . 

The other two obvious constraints read 

3 Q = 6^ + 2 c^a^ - Ub^ + 28 - 34 + + 259 + 2 fo^c^ > 

and 

-R = a^c^ - 10 b^(? + 30 c^a^ + 225 - 30 - - 25 6^ - 225 - 150 6^ > . 

The last constraint needed to define T>^^^ is then given by eq. ([8]). 

For illustrative purposes, the latter, purely algebraic description of the geometric 
shape of the boundary dT>^^^ may be complemented by the explicit evaluation of all 
the six energy levels (say, < Ex < ■ ■ ■ < E^). The easiest answer is obtained at 
the EEP singular couplings where the complete confluence of all the energies takes 
place, Eq = El = . . . = E5. 

A clear insight in the structure of the spectrum remains available in the small 
EEP vicinity as well. For example, the innermost pair of the energies E2 and E^ can 
coincide at dV'^^^ (and, subsequently, complexify out of T)^^^) while the remaining 
two doublets remain real, or vice versa. In the former scenario we encounter the 
pairwise coincidence of E2 = E3 = at dV^^\ In the latter case we may abbreviate 
Eq = El = —Az = —E4 = —E^ and compute the unknown quantity z G (0, 1). 
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Alternatively, the two outermost levels [viz., Eq = —E^ = —y^Fy where y G (0, 1)] 
can stay real while the confluence only involves the two internal energy doublets at 
the shared values of E^ = E^ = 2 ^/x = —Ei = —E2 where x G (0, 1). This is the 
most complicated example where one has the relation 

{s-Sraax) [s - Srainf = s^-{^2x^ + 25 y^) s^+(256x^ + 800 xV) s-6400xV = 
which deflnes a sub- surface of dV^^^ 



5.3 The pairwise confluences of the exceptional points 



Let us now return to the J = 3 option where the surface of the three-dimensional 
domain V^^^ can be visualized as composed, locally, of the two smooth sub-surfaces 
which intersect along a certain double exceptional point (DEP) curve. In terms of 
the single free parameter z, the DEP secular equation degenerates to the formula 
E^ {E + Azf {E-Azf = Q obtainable from eq. ©, 



(4.) 



~2{AzYs' + {Azfs = 0. 



(28) 



It is fortunate that the necessary analysis can still be performed non-numerically since 
equation fl28l) is easy to compare with the true secular equation ([5]) with coefficients 
given in paragraph 15. 2[ As long as the factorizable coefficient at s° must vanish, we 
get the first DEP constraint 



ac 



+ 15a+ (l5 + c^ + 56^ 
so that we may eliminate 



ac 



+ 15 a- (l5 + c^ + 5 6^ 
15 + c2 + 5 62 







c2 + 15 

In the quadrant of the a — h — c space with the positive a, the plus sign must be 
chosen, 

5 62 

a = 1 + 



c2 + 15 



Thus, we have 3 > a > 1 in the closed formula for 



62 = (c2 + 15) (a - 1) 
5 



(29) 



or, alternatively, for 



562 
a — 1 



15, 
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This result is to be complemented by the other two relations 



3Q{c,b,a) = 32z^ , 



R{c,b,a) = 128 2^ 



A straightforward elimination of z"^ 



gives the second DEP condition 



-66 - 36 6^ + 4 c^a^ - 189 + 252 0^-4 fo^a^ - = 



with the two compact roots 




2 - 33 - 2 6^ ± 2 Vc^ + 30 - 2 b'^c^ + 225 + 24 62 + . 



The acceptable one must be non-negative. For this would mean that 2 0^ > 
33 + 2 6^ while, at the same time, 63 + 12 6^ > 84 c^. These two conditions are 
manifestly incompatible so that we must accept the upper-sign root which is 
automatically positive for all the large 2 0^ > 33 + 2 6^. It also remains positive for 
all the smaller constrained by the requirement 



84 > 63 + 12 b^ . 

After the insertion of the definition fl29|) of 6^ we arrive at the final formula 



It defines a manifestly non-empty domain of parameters at which one encounters the 
pairwise confluences of the Kato's exceptional points. 

6 Summary and outlook 

There exist two reasons why our knowledge of the physical domains V is relevant. 
Firstly, their boundaries dT> are marking the breakdown of the reality and observ- 
ability of the spectra {En}- Secondly, these boundaries also represent a regime where 
the matrices H^^^ cease to be diagonalizable. Thus, it is the simultaneous degeneracy 
of the energies and of the wave functions which characterizes the physics near the 
boundary &D of the dynamical domain where the quantum system starts to become 
unstable. 




(30) 
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6.1 The onset of instabilities along dV 

In Landau's textbook ^22j on Quantum Mechanics the emergence of an instabihty 
of a quantum system is exemphfied by a particle in a strongly singular attractive 
potential V{x) = — G^/|xp. At the critical value G"^^^^^ = —1/4 of its strength one 
encounters a horizon beyond which the particle starts falling on the center. Vice 
versa, the system remains stable and physical on all the interval V = (—1/4, oo) of 
couplings G. From the pragmatic point of view the Landau's example is not too well 
selected since the falling particle should release, hypothetically, an infinite amount 
of the energy during its fall. A slightly better textbook example of the loss of the 
stability is provided, therefore, by the Dirac's electron which moves in a superstrong 
Coulomb potential. In the language of physics, particle-antiparticle pairs are created 
in the system beyond a critical charge {Z(^rnax) = 137 in suitable units j23]). 

What is shared by the above two Hermitian sample Hamiltonians is that they 
are well defined in a certain domain V of parameters while they lose sense and ap- 
plicability for parameter(s) beyond certain horizon. On a less intuitive level, similar 
situations have been studied by Kato [21]. He considered several finite-matrix toy 
Hamiltonians H{X). He paid attention to the unphysical, complex values of A and 
deduced that the related (in general, complex) spectra En{X) change smoothly with 
the variation of the parameter A unless one encounters certain critical, exceptional 
points A(^^). 

In the above context, certain carefully selected non- Hermitian examples seem 
to be able to offer the best illustrative examples in the physical stability analysis. 
In addition, our low-dimensional non-Hermitian tridiagonal matrices also path the 
way to a combination of mathematics and physics. In particular, in the context 
of pure mathematics, our present set of the solvable examples H^^^ with < 11 
enabled us to construct and study the real version of the Kato's exceptional points 
Xi^P) g dT>^^\ In parallel, the abundance of the J free parameters in our models 
looks more suitable for phenomenological purposes. In particular, we believe that a 
systematic characterization of a collapse of a realistic physical system could make use 
of these specific models offering an important link to the possible future classification 
and, perhaps, typology of quantum catastrophes [20] . 
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6.2 Degeneracies along the horizons dV 

The obvious theoretical appeal of the problem of stability may be perceived as one of 
the explanations of the recent growth of popularity of the PT— symmetric and, more 
generally, 77— pseudo-Hermitian Hamiltonians with promising relevance in quantum 
field theory [121 [23] and in quantum physics in general [HI [20] . Although the origin 
of the latter ideas can be traced back to the very early days of Quantum Theory 
|25j . the feasibility of its separate implementations have long been treated as a mere 
mathematical and/or physical curiosity (cf., e.g., [T3| [27] for illustration). 

One of the serious technical shortcomings of the PT— symmetric and other similar 
models is that their spectra are real (i.e., observable) in domains V with, sometimes, 
very complicated and strongly Hamiltonian-dependent shape of their EP boundaries 
dV. For an uninterrupted development of their study it may prove very fortunate 
that an explicit analytic description of the horizons dT> has been shown available 
here for all the matrix PT— symmetric chain models iJ^^^ with < 11. 

In this context it is particularly important that several recent microwave mea- 
surements [28] confirmed the observability of the abstract Kato's exceptional points 
in practice. These experiments re-attracted attention to the theoretical anal- 
yses of the EP horizons, say, in nuclear physics where many nuclei can, abruptly, 
lose their stability [TOl [29]. The growth of the role of the EPs may be also detected 
in the random-matrix ensembles with various interpretations [30] and in optical sys- 
tems (where EPs are called degeneracies [31] )• In classical magnetohydrodynamics 
the Kato's exceptional points may even happen to lie inside the domain of accept- 
able parameters, separating merely the different dynamical regimes of the so called 
a^— dynamos |5]. 

In all these contexts, our present completion of our recent studies of EPs may 
find its future role and relevance as a classification tool offering a deeper geometric 
understanding of the structure of the domains V (H). Basically, our results seem to 
indicate an efficiency of a combination of the methods of algebra (e.g., of solvable 
equations) and analysis (offering, e.g., the optimal parametrizations of elementary 
curves and (hyper) surfaces) with the computer- assisted symbolic manipulations and 
with perturbation expansions. Perhaps, our explicit verification of the complemen- 
tarity, compatibility and productivity of these methods could also lead, in a not 
too distant future, to the development of an explicit control of the stability of the 
systems, mediated by some purely algebraic tools of control of parameters in phe- 
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nomenological quantum Hamiltonians. 
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